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AN EQUATION THAT DESCRIBES MATERIAL OUTGASSING 
FOR CONTAMINATION MODELING 


Thomas M. Heslin 

Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

For many years, molecular contamination from material outgassants has been a continuing 
problem in the construction of spacecraft. Instruments with cooled detectors and optical 
instruments that operate in the ultraviolet light region are particularly susceptible. In many 
cases, contamination sources are at or near room temperature, while condensing surfaces 
are at a lower temperature. 

Prediction of the time- and temperature-dependence of room temperature outgassing rates 
is complicated because rates are usually so low as to make direct measurement impractical,* 
and also, the lack of molecular homogeneity in most materials complicates application of 
the Langmuir equation. The Langmuir equation relates the outgassing rate and the vapor 
pressure of a homogeneous monomolecular weight material during exposure to high vacuum. 
These two difficulties were encountered in a contamination problem on a series of very high 
resolution radiometers (VHRR’s) that were flown on weather satellites (ITOS D through G) 
for the National Oceanic and Atmospheric Administration (NOAA). Efforts to solve the 
contamination problem resulted in the derivation of an equation for a material’s vapor pres- 
sure that can be used in contamination modeling. Three examples of this type of application 
are given in the Appendix. 

DERIVATION 

Langmuir’s equation for the outgassing rate of a material exposed to high vacuum is 

dW r + 

— = 0.439 An/M/T P v J (I) 


*Direct measurements are impractical either because of the length of time involved or because the room temperature vapor 
pressures of outgassants are close to or less than the pressures obtained in outgassing test equipment. 

f Refer to the Appendix, Case 1, for application of Langmuir’s equation to a vacuum microbalance system and for the 
conditions under which use of equation 1 is justified. 


where 


dW 

dt 

A 

M 

T 

P v 


Weight loss rate, mg/sec 

Surface area, cm 2 , of sample 
Molecular weight, g/mole 
Absolute temperature, K 
Vapor pressure, Pascal, at T 


P y , according to the Clausius - Claperyron equation is 


P = P e - (E/R)(1/T -i/V 

V O 

where 

P o = Vapor pressure, Pascal, at T o 
E = Heat of vaporization, cal/gmole 
R = gas constant, 1.98 cal/gmole K 
Thus, 


dW 

dt 


= 0.439A n /M/T 


— (E/R) (1/T - 1/T q ) 


( 2 ) 


For an outgassing plastic or adhesive, E and P o are functions of molecular weight M, which, 
in turn, is a function of time and temperature. D. F. Wilcock (Reference 1 ) experimentally 
determined vapor pressure as a function of molecular weight and temperature for mono- 
molecular weight factions of methylpolysiloxanes. Generalizing from Wilcock’s work, a 
relation between molecular weight, temperature, and vapor pressure was presumed to be 


P 


V 


c M -c, M/RT 
Cj e 2 e 


where Cj , c 2 , c 3 are material constants. 


(3) 
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Substituting the right side of equation 3 for P y in equation 1 and defining three new constants 
y, E q , and M q such that e M (c 2 ~ c 3 /RT o ) = y M , c 3 = E o /M o , and c t = P o /y M o , results in the gen- 
eralization of equation 2, which is 


dw , — M - M -(E M/RMJ (1/T - 1 /T ) 

— = 0.439An/M/T P Q y ° e ° 


Extending this generalization to the Clausius-Claperyron equation gives 


p = P v M ' M ° e - (E o M/RM o> <‘/ T - »/T 0 ) 

V o ^ 


(4) 


(5) 


DISCUSSION 

M o is taken to be the molecular weight of the lightest faction of a given species and E o is its 
heat of vaporization. The factor P Q y M * M o gives the variation of vapor pressure with molec- 
ular weight at T o and the expression E o M/M q gives the variation of the heat of vaporization 
with molecular weight. The constant y must be less than 1 in order for vapor pressure to 
decrease with increasing molecular weight at T o . 

If a material were actually composed of a single generic species such as Wilcock’s methylpoly- 
siloxanes and its molecular weight distribution were given by 12, then there would be one 
expression such as equation 5 for each molecular weight, M, represented in S2. The total 
vapor pressure of the outgassants would be the sum of these individual partial pressures 
defined by equation 5 and the composition of. would determine how long the various 
partial pressures would exist. The problem becomes complicated at this point because the 
composition of S2 would have to be time-dependent. This is because some molecular weight 
factions would be exhausted before others as a function of the original amounts which were 
present and their rates of depletion. Also, if reaction or polymerization occur that change 
the distribution of molecular weights, these would have to be represented in €1 as a time- 
temperature function. 

In the outgassing of a real material, a number of generic species are usually evolved simul- 
taneously; therefore, the system becomes even more complex and makes the problem ex- 
tremely unwieldy analytically unless some simplifying assumptions are made. 

It has been found that an equivalent (from an outgassing point of view) material can be de- 
fined that allows the problem to be dealt with on a phenomenological basis. Using this 
approach a functional form of M is defined that makes equation 4 fit observed outgassing 
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data with the values of E and P equal to the constants for water at T . The most success- 

o o A o 

ful attempt at defining M in this manner has been with equations of the form 

M = M o X f(W/W o ), M o = 18 (6) 

where W q is conceptually equivalent to the “active mass” referred to in Papazin’s work 
(Reference 2) and f(W/W Q ), is a continuous function that increases with increasing fractional 
weight loss, W/W . The “active mass” W Q is an increasing function of increasing temperature. 
An equation of this type describes an equivalent material whose outgassing initially consists 
mainly of water. It assumes a material that predominantly outgasses vapor with a single molec- 
ular weight after a given weight loss at a particular temperature and it assumes that the 
molecular weight of the outgassed vapor continuously increases as sample weight decreases. 
Time is not explicit in equation 6 and there are no time-temperature, cross-product terms; 
therefore, the equation does not account for polymerization or reaction in time to form a 
different molecular weight distribution in the material. The extent to which the outgassing 
of a material cannot be described by an equivalent material such as described in equation 6 
is measured by the extent of the variation in the constants in the equation for f(W/W ) when 
it is fitted to observed outgassing data. 
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APPENDIX 


APPLICATION OF EQUATION 5 TO THREE COMMON 
OUTGASSING SYSTEMS 


Let x be the net weight of material lost from an outgassing sample per unit area. The rate 
of this weight loss is equal to the difference between the rate of evaporation, G (g/cm 2 sec) 
and the rate of condensation, C. 


dt 

The rate of evaporation from the material, according to the Langmuir equation (Reference 1 , 
p. 1 8), is 


G = 4.39 X 10‘ 4 \/m7t P y 


where P y is the vapor pressure (Pascal) of the material with molecular weight M (g/mole) 
being evaporated at temperature T (K). 

The rate (g/cm 2 sec) at which molecules of molec ular weight, M, strike the surface of the 
outgassing material is equal to 4.39 X 10' 4 y/M/T P; where P is the material’s partial pres- 
sure in the vicinity of the sample (Reference 1, p. 14). Assuming a sticking coefficient of 1, 
which is justified for most room temperature outgassants* (Reference 1, p. 18), then 


C = 4.39 X 10- 4 \/M/T P 


Thus 


— = 4.39 X 10' 4 n/M/T (P y - P) (A-l) 

Having established this, most outgassing contamination problems consist of determining P 
in terms of P y from constitutive equations based on the conservation of mass, and deter- 
mining P y experimentally. 


*The sticking coefficient referred to here is not the same as the sticking factor discussed in Reference 2 . The sticking 
coefficient referred to here is the probability that a molecule will strike the surface and stay there long enough to accomo- 
date thermally to that surface rather than rebound specularly from the surface. The sticking factor defined in Reference 
2 is a relative deposition rate normalized to the rate of 253 K (-20° C). 
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Case I. Condensing Surface in Series with a Pump 


This problem is outlined in the following flow diagram: 



The meaning of the flow diagram is that all the molecules outgassed from the sample must 
first condense on the room temperature collector and then be reemitted before being lost 
from the system. An example of this system is an Ainsworth balance often used to conduct 
outgassing measurements. Material outgassed from the sample goes to one of three places. 

It either returns to the sample and condenses as indicated by the P term in equation A-l , 
or it strikes the relatively cold chamber walls and condenses, or else it is reemitted from the 
walls and finds its way to the diffusion pump with its liquid nitrogen (LN 2 ) cold trap, and 
leaves the system. The latter two possibilities are expressed as follows: 


dx 

A 

dt 

material 
lost from 
sample 


4.39 X 10- 4 Aj n/M/T (P - P w ) + 

material condensed 
on room temperature 
chamber walls 


— - (P - PJ 

material lost to 
vacuum pump 


(A-2) 


where 

A = the sample area (cm 2 ) 

A, = the area where material condenses at room temperature (cm 2 ) 

P w = the vapor pressure (Pascal) of material with molecular weight M at the wall 
temperature 

P b = the vapor pressure of material with molecular weight M, in the diffusion pump 
cold trap 

F = the conductance (liter/sec) between the point where vapor pressure is P w and 
where it is P. 

b 

R = Gas Constant = 8.294 X 10 3 Pascal - liter/mole K 
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Solving equation A-2 for P and substituting into equation A-l and solving for dx/dt gives: 


dx 

dt~ 


4.39 X 10‘ 4 A lV /M/T 
(Aj + A) 


MF 

(P - P ) + 

1 v w RT (A + A) 


(P w " P b> 


(A-3) 


Equation A-3 expresses the same thing that equation A-2 does and interpretation of the two 
terms involved is the same except now the equation is expressed in terms of the vapor pres- 
sure of the material. 

Equation A-3 can be applied to an Ainsworth balance system used at the Goddard Space 
Flight Center (GSFC). A schematic of the system is shown in figure A-l. 

Material is observed to condense on the inside of the quartz chamber wall along a 5- to 1 0-cm 
long area immediately adjacent to the heater. Thus, if the length of a condensation area 
is taken as 10 cm 


^ =2 X 2n (10) = 479 cm 2 

The conductance between a point just outside the furnace and the diffusion pump is complex. 
Davis (Reference 3) has calculated F for cylindrical elbows and other geometries that prevent 
line-of-sight transmission. However, this is not needed because the mass flow for material 
that did not condense at room temperature just outside the furnace is the same at any point 
along the line to the diffusion pump, or else an accumulation would be observed. Therefore 
its value can be determined at any point along the line. A measure of this flow is given by 
the expression: 


MF X ( P w " V 
RT (Aj + A) 


(A -4) 


where P b is the material’s vapor pressure at some point between the area where its vapor 
pressure is P w and the diffusion pump cold trap where its vapor pressure is P b . The constant 
Fj is the conductance in the area between the points where the vapor pressure is P w and 
where it is P b . Expression A- 4 can therefore be substituted for the second term in equation 
A-3. A pressure measurement is determined at the base of the first elbow in the system by 
an ionization gage that is calibrated for air. It typically reads in the 1.33 X 10' 4 Pascal range 
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15.24 CM DIA. 


25.40 CM LONG 
10.16 CM DIA. 


33.02 CM LONG 
10.16 CM DIA. 


22.86 CM LONG 
7.62 CM DIA. 


QUARTZ TUBE 
JUST ABOVE 
HEATER 


Figure A-1 . Schematic diagram of an Ainsworth balance. 




soon after evacuation has started (figure A-l). The conductance from just above the furnace 
to this point is relatively simple and is given by (Reference 1, p. 88) 


d 3 /T\* 

Fj = 30.48 — y'jJJ ) liter/sec 


where d = the diameter of the tube (cm) and L is its length in cm. For a material with molec- 
ular weight 1 8 at 300 K in the Ainsworth balance 


Fj = 30.48 


(7.62) 3 /300\* 

(22.86) V 18 / 


= 2408 liter /sec 


The sample area is typically 20 cm 2 . Substituting expression A-4 in equation A-3 and using 
6.65 X 10‘ 4 Pascal for P b gives 


dx 

dt~ 


1.05 X 10' 4 (P y - P w ) + 3.49 X 10' 5 (P w - 6.65 X 10' 4 ) 


(A-5) 


for T = 300 K. 

Ordinarily A « A t and at temperatures greater than room temperature P y » P w and 
6.65 X 10' 4 Pascal, therefore equation A-3 reduces to 

= 4.39 X 10' 4 v/m/T P y g/cm 2 sec or ~ 1 .05 X 10' 4 P y (M = 18 g/mole, T > 300 K) 

which is subtan tively equation 1 in the main text. 

However, at room temperature P y - P w and equation A-5 points up a caution that must be 
exercised in applying equation 1 once the equation for P y is determined. Equation 1 , which 
adequately describes high temperature outgassing rates, cannot be used to give room tempera- 
ture rates for this type of system. As room temperature (300 K) is approached, the out- 
gassing rate for M = 1 8 is approximated more closely by 

Hy MF (P - P. ) . 

= ? — — 1- ~ 3.49 X 10‘ 5 (P - 6.65 X 10‘ 4 ) g/cm 2 sec (A-6) 

dt RT (Aj + A) 

Equation A-6 also indicates that materials with a room temperature vapor pressure less than 
P b ~ 6.65 X 10‘ 4 Pascal do not appear to lose weight in this system and can even appear to 
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gain weight slightly once very volatile materials have been removed and backstreaming pump 
oil starts to condense on a material that has a large surface area. This phenomeon has been 
observed with aluminum anodize, for example. 

For Case I systems, the room temperature collection term is generally of the most interest since 
in practical applications, other than the Ainsworth balance system, the room temperature 
collector surface usually represents an optic of some sort. If S represents this term then 


S = 4.39 X KT 4 Aj/M/T (P - P w ) 


(A-7) 


Solving equation A-l for P and substituting into equation A-7 and substituting equation A-3 
for dx/dt gives, on rationalizing 


S = 


4.39 X 10' 4 AAjTm/T 
(Aj + A) 


(P - P ) - 

V V V/' 


AjMF 


(P 

(Aj + A) w 


PJ 


Case II. Condensing Surface in Parallel with Pump 

This problem is illustrated in the following flow diagram: 



In this case the constitutive equation resulting from the conservation of mass is 


UA . / ivir 

A — = 4.39 X 10' 4 Aj yM/T (P - P w ) + — (P - P„) (A-8) 

at K1 


Solving equation A-8 for P and substituting into equation A-l and then solving for dx/dt gives: 

0% 


rl y 

= 4.39 X 10' 4 7m/T 

dt v 


(P y - P w ) A t RT Vm/T + (P v - P t 


MF“ 


-4^X-tb' 4 TdF + (A x + A) RT/M/T 
The room temperature collection term, S, is 

|(P V - P w ) ART n/m/T - (P w - P b )Hr39-X-WrnfF 


S = 4.39 X KT 4 A, /m/T 


^^Cib^F + (Aj + A) rt/m/t 

3;na \-np 


(A-9) 
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Case III. No Return Flux due to Container 


In this case there is no interaction of the evolved gas with anything but itself. To approximate 
this situation the outgassant is imagined to evolve into an expanding volume, V, whose radius, 
r, is given by the mean free path. The constitutive equation is therefore 


P 


nRT 

V 


Ax RT 
M V ’ 



r = L o + vt 


(A-10) 


where 


= initial mean free path 


= mean molecular velocity = 



Vz 


Substituting equation A-10 into equation A-l gives 


dx 

dt 


4.39 X l(r 4 yM/T 



x3ART 


4rrM (L q + vt) 3 


In this case if a small room temperature collector is introduced which returns relatively 
little material to the sample, then the room temperature collection term, S, is simply the 
product of the outgassing rate times a view factor assuming a sticking coefficient of one. 
The view factor is determined in the same way as for a radiative heat transfer problem. 
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